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Abstract 

Based on the Unruh effect, we calculate the critical acceleration of the Bose-Einstein 
condensation in a free complex scalar field at finite density in the Rindler space. Our model 
corresponds to an ideal gas performing constantly accelerating motion in a Minkowski 
space-time at zero-temperature, where the gas is composed of the complex scalar particles 
and it can be thought to be in a thermal-bath with the Unruh temperature. In the 
accelerating frame, the model will be in the Bose-Einstein condensation state at low 
acceleration, on the other hand there will be no condensation at high acceleration by the 
thermal excitation brought into by the Unruh effect. Our critical acceleration is the one 
at which the Bose-Einstein condensation begins to appear in the accelerating frame when 
we decrease the acceleration gradually. To carry out the calculation, we assume that the 
critical acceleration is much larger than the mass of the particle. 


1 Introduction 


The Bose-Einstein condensation is getting a lot of attention recently as a quantum fluid for 
the test of the analogy between sound waves in quantum fluids and scalar field fluctuation 
in curved space-times [Tj . Thanks to this analogy, we can expect new progress and insight 
which are difficult only in the gravitational analysis in the quantum gravity and cosmology. 
Further, the experiments of the quantum gravity and cosmology are difficult due to the 
required energy level and the scale of phenomena. However, the experiments of the 
condensed matters will be free from such problems to some extent. It is thought for these 
reasons that pseudo experiments of the gravity are possible in the quantum fluids through 
the analogy. 

In order to understand the quantum phenomena in the gravity and cosmology such 
as the Hawking radiation [2] and particle creation EQBI, the Unruh effect 0 El HP is 
important in terms of the role that the event-horizons play. 

The Unruh effect is a prediction that one moving in the Minkowski space-time with a 
linear constant acceleration experiences the space-time as a thermal-bath with the Unruh 
temperature, T v = ha/ (/I-ncks) ~4x 1CU 23 a/(cm/s 2 ) [K], where a is the acceleration. 

Now, various experimental attempts in the condensed matters to observe the grav¬ 
itational phenomena are being invented (see Ref.[9] for example). Particularly as for 
the experiments to detect the Unruh effect, there are attempts in Bose-Einstein conden¬ 
sates [11], graphenes [12] and Berry phases [[[3]. For other attempts see Ref. [IT], for 
example, and related references. 

We also address the issue of the Unruh effect in the Bose-Einstein condensation. 
Whether the Bose-Einstein condensation occurs or not is determined by temperature. 
We assume in this paper that the Unruh temperature exists in the constantly accelerat¬ 
ing system according to the Unruh effect mentioned above. At this time we can think that 
the Unruh effect affects the Bose-Einstein condensation. Although an enormous number 
of studies have been done on the Unruh effect and the Bose-Einstein condensation so far, 
these are performed separately and little is known about the Unruh effect in the Bose- 
Einstein condensation at this moment. Since both the Unruh effect and the Bose-Einstein 
condensation are very important in the fundamental physics, new understandings could 
be expected by combining the Unruh effect and the Bose-Einstein condensation. In this 
paper, from such a background, we calculate the critical acceleration for the Bose-Einstein 
condensation based on the thermal excitation brought into by the Unruh effect. 

Let us here explain the Bose-Einstein condensation briefly (for more details, see 
Ref. [TO] for example). The Bose-Einstein condensation state is the situation that all 
the particle stay in the least energy state uniformly owing to the Bose-Einstein statistics, 
and it appears as a phase transition that all the particles uniformly drop to the least 
energy state as the entropically-favored state due to the Bose-Einstein statistics at some 
time. (In the Bose-Einstein condensation state, since all the particles are in a state such 
as the least energy state, de Broglie wave of each particle becomes longer and even, and 
eventually the system itself becomes a de Broglie wave.) The situation that the effect of 
the Bose-Einstein statistics is dominant can be considered as the low temperature region. 
Actually it is considered at the low temperature region that the particle’s thermal motion 
energy is extremely small so that it does not excite the state of the particles from the 


1 


least energy state. Hence, the Bose-Einstein condensation state emerges in the system 
composed of bosonic particle at extremely low temperature. 

We calculate in this paper the critical acceleration at which the Bose-Einstein con¬ 
densation begins to appear in the accelerating system by the Unruh effect when lowering 
the acceleration. We mention our critical acceleration more precisely. We first consider 
an ideal gas in the Minkowski space-time at zero-temperature, where the gas is composed 
of particles described by a free complex scalar held at finite density. As the gas is now 
in zero-temperature, it can be considered to be in the Bose-Einstein condensation state. 
We then start to accelerate such an ideal gas uniformly. At this time, according to the 
Unruh effect, since each constantly accelerating particle composing the gas will experi¬ 
ence the temperature Tu = Ha/(2nckB) (a is the acceleration) in the accelerating frame, 
the ideal gas performing uniformly accelerating motion will experience the temperature 
Tu as a whole in the accelerating frame. Here, what the gas experiences thermal means 
that the space which the gas observes is filled with some medium that are performing 
thermal fluctuation. In this paper, we consider that there is no interaction between the 
medium and the particles composing the gas without the thermal excitation. Hence it 
is considered that, as growing the acceleration gradually from lower acceleration, in the 
accelerating frame, the dissolution of the Bose-Einstein condensation is observed eventu¬ 
ally, and finally the Bose-Einstein condensation disappears entirely at some acceleration. 
The critical acceleration we calculate is the acceleration at that time. 

We here mention the emergence mechanism of the Bose-Einstein condensation in our 
model and its critical moment. It is composed of three steps. We first calculate the 
effective potential of the held meaning the particles composing the gas, and then obtain 
the particle density by performing a derivative to it with regard to the particle’s chemical 
potential. Then we find that, if we decrease the acceleration with fixing the particle 
density to constant, either the particle’s chemical potential or the absolute value of the 
zero-mode of the held has to grow. Here, the acceleration plays the role of the temperature 
in the accelerating frame as mentioned in the above paragraph, and the zero-mode of the 
held can be considered to correspond to the least energy state. Hence the absolute value 
of the zero-mode is considered as the expectation value of the Bose-Einstein condensation 
state, and whether it is zero or not corresponds to its disappearance or appearance, 
respectively. We consider to start with the high acceleration situation where there is 
no Bose-Einstein condensation in the accelerating frame, and decrease the acceleration 
gradually. At this time, we keep the particle density to constant. Namely, the number 
of particle is always constant in our model. Then, corresponding to no Bose-Einstein 
condensation, the absolute value of the zero-mode should be zero while the we decrease 
acceleration. Thus the chemical potential must grow to keep the particle density to 
constant (step 1). However we find that there is an upper bound for the value of the 
chemical potential, which is need to avoid a divergence of the probability amplitude. 
Hence, when the chemical potential reaches the upper bound (step 2), there is no way 
but the absolute value of the zero-mode starts to grow (step 3), if we further decrease 
acceleration with keeping the particle density to constant. Thus, the moment that the 
chemical potential reaches the upper bound is the critical moment that the Bose-Einstein 
condensation appears, and this is the Bose-Einstein condensation in our model. 
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The situation where our analysis is actually performed is just before the Bose-Einstein 
condensation state starts to occur. To be more specific, our calculation to obtain the 
critical acceleration is performed in the situation that the absolute value of the zero-mode 
and the chemical potential are put to zero and the upper bound value respectively in 
the equation of the particle density as the critical moment mentioned above. Then the 
acceleration obtained from that equation is the critical acceleration. 

In our analysis, due to some technical difficulty mentioned later, we consider the 
situation m/a c <C 1 (m and a c are the mass of the particle and the critical acceleration, 
respectively), and take the leading contribution of this in our calculation, where m is 
the mass of the complex scalar particle considered in this paper and a c is the critical 
acceleration. As the acceleration is proportional to the temperature in the relation given 
by the Unruh effect, this can be read as ksTj 2 me 2 (T c is the critical temperature 
fixed with a c ), which is the relativistic situation that the particle’s thermal motion energy 
is much higher than its static energy at the critical moment. 

The result in this paper is hence applicable to the system of the complex scalar field 
with the critical temperature for the Bose-Eiusteiu condensation that is much higher than 
the mass. Since the kind of the field considered in this paper is a complex scalar field 
composing an ideal gas, the actual particle that can correspond to this paper is some ideal 
complex scalar particle system whose critical temperature of the Bose-Einsteiu condensa¬ 
tion locates in the region ksT/2 me 2 . 

As stated above, in this study we perform the calculation for the critical acceleration 
in the Bose-Einstein condensation. So far several kinds of the critical accelerations for 
the spontaneous symmetry breaking induced by the Unruh effect have been carried out. 

The critical acceleration for the chiral symmetry restoration was studied in Nambu- 
Jona-Lasinio model at zero and finite chemical potentials for quarks in refs. [15] and [16] 
respectively. In Ref. B33 the critical acceleration for the restoration of the spontaneous 
symmetry breaking of the Z 2 symmetry in the real scalar field theory was studied. 

There are also papers concluding that the Unruh effect does not contribute to the 
restoration of the spontaneous symmetry breaking [18]. The system in these papers is not 
the finite density, and their discussion is based on the effective potential. On the other 
hand, our Bose-Einstein condensation occurs from the relation between the chemical po¬ 
tential and the acceleration playing the role of the temperature. This is a different point 
between us and him. Because of this, their conclusion would not be applied to our study 
readily. We should take this issue to a future work. 

Lastly, we introduce interesting papers that will have some connection to this study. 
Ref. [19] concludes that a larger acceleration should enhance a condensate as compared to 
those in a non-accelerated vacuum. In Ref. [SD], although the background space-time is 
not the Rindler space, whether the Bose-Einstein condensation can occur or not is shown 
in an ideal boson gas model with a point-like impurity at finite temperature in some 
uniform gravitational force in each of D = 1, 2 and 3. Next, it is shown in Ref. ED that 
free massless scalar particles are detected by a constantly accelerating detector with not 
the Bose-Einstein distribution but the Fermi-Dirac distribution in odd dimensions, and 
this problem is solved in Ref. [S2]. 
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2 The model 


The model in this paper is a free complex scalar model at the finite density correspond¬ 
ing to an ideal gas in the constantly accelerating system with the acceleration a. The 
Lagrangian density is given as 

C = h 2 g' w d^cj)* d u (f> — c 2 m 2 (j)*(j) (1) 

with (j> = ( 0 i + i(f> 2 ) that means the field of the particle composing a gas, and p, u are 

the coordinate in the Rindler space explained below. 

In the Minkowski space-time as an inertial system, the background space-time of a 
constantly accelerating system is given by the Rindler space. The Rindler coordinate in 
our notation is given as follows: 


(Vi Pi Vi z ) = (V: Pi x -l) with x± = (y,z). 

This relates with the Minkowski coordinate (t,x,x±) as 

(t,x) = — (sinh—, ccosh—) = p(sinhp, c cosh 77 ), 
a c c 


( 2 ) 


(3) 


where the accelerating direction has been thought to be in the ^-direction. The Rindler 
metric in our notation is 


ds 2 = ( cp) 2 dr) 2 — d{cp) 2 — dx (4) 
In what follows, we use the unit system: c = h = ks = 1. 

The constantly accelerating one in the Minkowski space-time corresponds to the one 
moving along a line on a constant p in the Rindler space. The relation between p and 
a are p = 1/a as can be seen from eq.((3]). Since constantly accelerating one experiences 
the system as a thermal-bath with the Unruh temperature Tu (the acceleration are in the 
relation: T v = a/2n) by the Unruh effect, one moving along a line on a constant p in the 
Rindler space gets the temperature 


T v = l/2irp. (5) 

Hence, the gas in our study is considered to be in such a thermal-bath with Unruh 
temperature Tu in the accelerating frame. Here, although our gas is considered to be 
in a thermal-bath with some medium performing thermal fluctuation giving the Unruh 
temperature Tu in the accelerating frame, we assume that there is no interaction between 
our gas and the medium other than the thermal excitation. 

We can see from eq.(l5l) that varying p means varying the temperature. For this 
reason, how to interpret the results would be unclear if the four-dimensional space-time 
integration including the p-integration was performed. 

Let us here turn to how this problem is handled in other papers on the critical ac¬ 
celeration given by the Unruh effect. In R.efs. [T5l fl 6 l ITT] , the action is given with the 
four-dimensional space-time integration including the p-integration. This point is the 
problem. However, their calculation is once performed using the Green’s function that 
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has dependence on the p-direction, the coordinate p is treated as a constant in the cal¬ 
culation of the effective potential. As a result, the p-integration becomes just a volume 
factor in the calculation of the effective potential. 

The difficulty in the treatment of the coordinate p is also mentioned at the chapter of 
conclusions and discussions in Ref. [23j in the context of the study of the Larmor radiation 
with the correction rooted in the Unruh effect. 

In our analysis, as well as Refs. [Tol IIBI RTj . we obtain the Green’s function that has 
dependence on the p-direction. However, the space-time integration in our action is given 
without the p-integration. As a result, p is a parameter in our model. 

We can see at this time that the action is needed to be multiplied by a quantity with 
the dimension of length so that the action becomes dimensionless. To this purpose, we put 
dp in our action. This means that the process to obtain the path-integral representation 
given in eq.([B]) for the probability amplitude from the operator formalism representation 
has been performed with fixing p. If the analysis in Refs. |T51 IBj i lTl is performed again 
by our way, there is no difference in the final result. 

There are four regions separated by the event-horizons in the Rindler space. The 
region treated in this paper is the right wedge only. 


3 The effective potential 

3.1 Performance of the path-integral 

We start with the probability amplitude: 


Z = / T > 7 t ?? 'E , 0 exp 


i / d 3 x dp'y ( + n 71 *d^cp* — {H — pq) 


( 6 ) 


where 7 = yf —det g Mi , and d 3 x = dr]d 2 x_i_, and p and q are the chemical potential and 
the density of particles, respectively. For the reason mentioned in Section.2, to fix the 
acceleration and Unruh temperature, the integration in the p-direction is not included in 
the space-time integration in the action. As a result, p is a parameter in our analysis, 
and the Unruh temperature that the gas experiences is given according to the relation in 
eq. At this time, we can see that the action is needed to be multiplied by a quantity 
with the dimension of length so that it becomes dimensionless. We put dp to this purpose. 
However the space-time integration is the one without the p-direction. Explanation for 
7 r v and H is in what follows. 

(tt" . tt'") = (tTTt- ajpj) are the momenta given as 

P\^) = (g^d^\g^d^). (7) 

V. is the Hamiltonian density given as 

H = n v d v (j) + 7i v * — £ 

= n 71 *^ — g^di^djcf) + m 2 4>*(j), ( 8 ) 


5 







where is defined as 7 ( 7 ^ + 2712 ). Correspondingly, the functional integration 
measure of tt v changes as = Vn^V^ri- From eq.([7|), we can see 

fa, 2 = g™d v fa, 2 . (9) 

It turns out that the conserved current associated with the U(l) global symmetry in 
our model is obtained as 


J» = -ig^fad.fa -fadfa). ( 10 ) 

Using this J M , the integral of conserved charge density can be written as 


J d 3 x dp'yq — j d 3 x dp^ J v = j d 3 x dp 7 ( — 7t 2 fa + fa fa) . 
Substituting eqs.(JH]) and (fTTjl into eq. (J 6 ]) , we can obtain the following Z : 

Z = j 'Dni n 'Dii2' n T > (i) exp - j d 3 x dp 7 j 

^{fa^iv - 2(d v fa + ufa) fa) - ^(fan 2r , - 2 (d v fa + P fa) fa) 


2 

^(( difa ) 2 + ( difa) 2 - m 2 (</> 2 + <j>. 


where i, j = x±. Here, in the above, we perform the following rewritings: 

fan lv - 2 (d v fa + p fa) fa = g m |(tti^ - (d v fa + pfa)) 2 - (d v fa + p 0 2 ) 2 }, 
fan 2v - 2 (d v fa + pfa) fa = c/ w |( 7 t 2j? - (d v fa + pfa)Y ~ {fa fa + pfa) 
Furthermore, we redefine the fields as 

TTl r) ~ (dfal + Pfa) ~t 7 Ti, ? , 

^2r) — farjfa + p fa) —> 7T 2t? . 

As a result, we can write Z as 


Z = C Vcftex p 


L 2 


d 3 x dp 7 [g^finfa + pfa) + (difa) 2 

+ g m (fafa + h0i) 2 + (difa) 2 — m 2 ((pi~ + 


( 11 ) 


( 12 ) 

(13a) 

(13b) 


(14a) 

(14b) 


(15) 


with C= T>'Ki v 'Dn 2 r, exp 


. Performing the path-integral 


i J d 3 x dp^g riri {('Ki 71 ) 2 + ('K 2 r 1 ) 2 ) 

of 7 r lj? and 7T 2r? formally, we think that C become some factor. We ignore C in what follows. 
As a result, with some straight forward calculation, we can write Z in the following form, 


Z = / F>0exp 


d 3 x dp~i (faGfa + faGfa + 2 g m p (fadfa 1 - fad v fa) 


(16) 
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where G = d 2 + g lj d^dj) + m 2 — g^g 2 . 

Let us now rewrite the real and imaginary parts of the held into a convenient expression 
for the analysis of the Bose-Einstein condensation. As we have written in the introduction, 
the Bose-Einstein condensation can be considered as the situation that all particles get 
to be the least energy state. The least energy state can be considered as the zero-mode of 
the held in our model, and the situation that all the particles are in the least energy state 
can be considered as the condensation of the zero-mode. Hence, a convenient expression 
for the analysis of the Bose-Einstein condensation in our analysis is the one in which the 
zero-mode is separated as 

0i = y/2 a cos 6 + 0i, (17a) 

02 = y/2a sin# + 02 , (17b) 

where a plays the role of the expectation value of the condensation of the zero-mode, 6 
is a phase in the zero-mode and 0i, and 02 are the non-zero modes. Depending on just 
before or after the condensation starts, a behaves as follows: 


a = 0 : before the condensation 
a^O : after the condensation 


At this time, 0iG 0i and 0 2 G 02 can be written as 

0iG 0i = 2 a 2 {m 2 — g vv g 2 ) cos 2 9 + 0 X G 0i, 
0 2 G 02 = 2 a 2 (m 2 — g vri g 2 ) sin 2 6 + 0 2 G 0 2 . 


As a result, we can rewrite Z into 


Z = 


exp 


— ia 2 d 3 x dp 7 (m 2 — g^g 2 ) 


T >0 ex p 


l -Jd 3 xdp t(A 4)( 2 s ij, Sjj 2,J '' G lldv 

V2 a cos 9 f G (0i + 0 2 ) + (0i + 02 ) G 


(18) 

(19) 


( 20 ) 


Our analysis will be performed just before the condensation as mentioned later. For 
this reason, we put a = 0 in what follows. As a result, we can write Z as 


Z = / D0 exp 


1 

2 


03 , (1 z G - 2 g^gdr, 

dxdp 7(0i 02 )( 2gri>d Q 


( 21 ) 


Then, we can see from the form of G given below eq. fllbl) that there should be the 
condition: 


m 2 - g v y 2 > 0. (22) 

Otherwise the path-integrals of Z diverges at the configuration that all the momenta are 
zero. Further, the above relation gives the upper limit of the chemical potential for given 
a mass and an Unruh temperature. 
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Performing the diagonalization as 


Z = 


I V([> exp 
Vp exp 


-J<Pxdp-,{h fo)UU-'^ 2grl>a 


a -2 ^a, ]uirl 

y rj 


01 


o 

G - 2 gWfj, d v 


( 23 ) 


where U = f * * j is a unitary matrix defined to perform the above diagonaliza¬ 
tion, and correspondingly | = U~ 3 ( - 1 |. At this transformation, the functional 

V02 / V 02 / 

measure is also transformed, which we have described as Deft —> V'cp'. However, since U is 
a constant unitary matrix, the difference between X >0 and contribute only to some 
constant factor in the path-integral, and we ignore it in what follows. 

Performing the path-integral, we can obtain 


Z = Det [{G + 2g™p d v ) (G - 2 g^p d v ) 


- 1/2 


(24) 


In the above, with regard to the treatment of dp, we consider that the integration of this 
has been performed by assigning a value at one point. Hence, the effective action W 
defined as Z = exp iW can be written as 


W = l - Log Det ((G + 2 g^p d r] ) (G - 2> d v )) 

= \ Tr Log ( (G + 2 \g™p d v ) (G - 2 g™p d v ) ) 

= JWr ( Log ( a ' + ^ + m2 + 

+ Log(<9j + 'y~ 1 g l3 di{yd j ) + M 2 - 2ig m p d v )^j, (25) 


where M 2 = m 2 — g vv p 2 and V = J d 3 x 7 is the volume for the (r),x±) space-time, which 
appears from the rewriting of a functional trace into an integration: 


Tr —>• 



(26) 


with (^) 3 = f dk 3 and V = L 3 = f d 3 x 7 , where L means just length in each space 

for the rj, ^-directions, and yj —det g a b = \J—g PP det g a b = 7 (a,b — rj, x± except for p, 
and g pp = 1 ). d 3 k = dujd 2 k± with k± = (. k y , k z ), which are momenta corresponding to the 
coordinate (rj,x±) as in eq. (132]) . 










Defining 5± = g m d^ + 7 1 g lj ^(7 ct,) ± 2 ig V11 pd 11 , we further rewrite IF as 

W=t 2 V Jj^( L ° g ^ ++ M + L ° g (' 5 - + M2) ) 

= ) h/FT ( / dA 2 (5 + + A 2 ) -1 + Log 5+ + ^ d.A 2 (5- + A 2 ) -1 + Log . 

(27) 


Here, we can ignore f t ||?3 Log<5± for the reason in what follows: First, since we will 
perform the derivative with regard to the chemical potential to obtain the particle density 
at the end, we look at only the part concerning the chemical potential as 


f dk 3 

(2 

r dk 3 

(2vr) 3 

* dk ' 


Log5± 


Log (d 2 + 7 1 g t 3 d i ( 7<9/)) + / dk r 


dk\ 

( 2 ^ 


Log 1 ± 


2 ig^iidr, 


<% + 'y~ 1 g ii di('yd j 


Log ( 1 ± 


2 ig^fJLdr, 


(2vr) 3 V <% + 'Y~ 1 9 ij di('yd j ) 


(28) 


In the above, we have described the integral f dk 3 separately as J dk v f dk\ . Then, 
replacing d t with i ki , from the fact that J dk k log ( 1 + c ^. 2 ) = 0 (ci j2 are some constants), 
we can see that the //-dependent part vanishes in the /^-integrals. 

Finally, we can write the effective action W in eq . 112711 as 

H, = ^V(i 5 P A2(S+ + S - ) ’ (29) 


where D± is 


7^(4 +A 2 ) 1 = D ± = D ± (k VlP ,k ± ). (30) 

In the above, we wrote the arguments concerning the momenta as k v and kj_ despite 
that these are given by the differential operators in the actual expression, which may be 
allowed. 

3.2 Calculation of D± 

From the definition in eq. fj30|l . we can have the following equation: 

5 A (x — x') =7 (<5± + A 2 ) D±{x — x') 

=7 {g m d^ + 7~y j <9i(7<9j) ± 2/ g m gd v + A 2 ) D±{x - x') (31) 

with D-i(x — x') dehned as 

D±{x ~x')= jD± e* («(W)-fcL(*L-^)). ( 32 ) 
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We can rewrite eq. (l3Tj) as 


S(p - p') 
f d 3 k 


dk J (u (r)-r)')-k±(x±-a! ± )) 


(2tt) 


(2tt) 3 

d 3 k 

( 2^)3 


7 + T'gVSihdj) ±2ig m p8 n + A 2 ) D ± e‘ 

7 (y 2 ^ - yyy,) - (6g + d;) ± 2 ig-va, + a 2 ) 5 ±e ‘ 

(33) 


where = dujdk\. As a result, we can obtain the equation which D± should satisfy as 

7 ~ l d{p ~ p') = ( - P~ 2 u 2 ~ p~ X d P ~ d 2 p + k 2 y + k\ ± 2p~ 2 pu + A 2 )D±. (34) 

Finally, we can obtain the equation that determines D± as 

(j^ + f ll)D± = -pS(p-p'), (35) 

where T = p 2 d 2 p + pd p - p 2 K 2 with k 2 = k 2 + A 2 , (fc 2 = A : 2 + k 2 z ), and Q 2 ± = uj(oj t 2/i). 

Let us now determine _D± from eq. (1351) . First, we can see that T can satisfy the 
following differential equation: 


with 


TG x (p, k ) = (A) 2 0 A (p, k). 


©a (p, k ) = CxKix(np), Cx = - \/2A sinh(7rA), 

7r 


(36) 


(37) 


where the values of A are positive real numbers, Kix(np) is the second kind modified 
Bessel function, and C\ was obtained from the equation that the second kind modified 
Bessel functions satisfy, 


f °° dp 

/o P 


©a '(/£>, k) @x(p, k) = 5{ A'-A). 


(38) 


We then assume that D± can be written by taking © A (p, k) as the bases for the 
orthogonal directions labeled by A as 


D + = 


d\f\,±(ui,p')Q x (p,k), 


(39) 


where the above means an orthogonal base expansion of Q x (p,k) labeled by A, and 
fx,±(w,p') are the coefficients of each direction. If we can find fx,±(w,p') that satisfies 
eq. (|35]) . this assumption is right. We now obtain such fx,±(w,p'). 

We can see that now we can have the following two independent equations as 


dp 

P 

dp 


D± T ©a(p, k) = - A 2 


dp 

P 


D± ©a (p, fc), 


dp 


— ©a (p, k ) TD ± = -O 2 / 0 A ( p , k ) D± _ © A ( p ', k ). 


P 


P 


(40a) 

(40b) 


10 






Here, in obtaining eqs. (l40ap and (14 Ob ll . we have used eqs. (l36|) and (T35T) . respectively. 

Eq. (l40all and eq. MObp are equivalent to each other. Actually, this equivalence can be 
seen easily by what eq. (l40ap and eq. (|40bjl can be represented as (D±\ .F|© A (p, k)} and 
(0 A (p, k)\F\D±), respectively. Then, by subtracting these two equations each other, 

0 = (- A 2 + fi|) / C ^D ± 0 A (p, k) + 0 A (p', k) 

J P 

= (-A 2 + fi|) f ^- Jd\'h'^u,p')G X '(p,k)-G x (p,k) + G x (p',k) 

= (-a 2 + h|) Jdx '• J ^e X '(p,k)e x (p,k)- fx, ± {u:,p') + ® x {p',k) 

= (- A 2 + fi|) f d\' ■ 6( A' - A) • / A ',±(w, p') + 0a(p', k) 

— (— A 2 + p' ) + 0 A (p / , k). (41) 


In the above, we have used eqs. (1381) and (l39p . From the above, f\.±{oj,p') can be deter¬ 
mined as 


/a,±(w,p') 


BaG o',k) 
a 2 -ni 


e A 77) 

A 2 - fi 2 ■ 


(42) 


In the above, we have regarded Ol and fi 2 as the same each other, and written as 
fi 2 = fi 2 = fi 2 . Because the integrations of fx^+(uj,p') and f\_((jj,p') with regard to oj 
in eq. (1T21) do not produce any difference under the transformation u —> —uj. Namely, 

n ( T — t') — I — < ^ ± [p\ Ba(P ; k) 0 A (P, k) j (u (r)-T)')-k±(xj_~a! ± )') 

+ [X X) ~ ] (2vr)3 'J dA A 2 -(u(o;-2 p) 6 

[dud k± f .. Oa(p j kdj O A (p, k') ; (u> (ri—r]')—k±(x±—a! ± )'j 

J ( 27 t ) 3 J X 2 -U}(cj + 2p) 

= D_(x — x'). ( 43 ) 


As a result, we can put = fi 2 . Correspondingly, there is no difference between D + 
and £)_. Hence, we put D± = D and D± = D. 

Now, substituting eq. (1421) into eq. (j39p . we can write D as 


D 



Q A (p / ,fc)0 A (p, k) 
A 2 - Q 2 


(44) 


Finally, the effective action W can be written as 

r pu3 rM 2 _ 

w ^Sw?L dA2 °- 


(45) 
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3.3 Calculation with the Euclidianization 

Now that D± = D have been obtained, backing to eq. (l29p . we perform the integrations. 
We start with performing the Wick rotation toward the 77 -direction. At this time, the 
metric given in eq.flTJ) changes to 

ds E = p 2 drf + dp 2 + dz? ± . (46) 

We can see from this form that the 77 -direction is 5' 1 -compactihed with the period /? = 2tt. 
Then, the following replacement arises: 


with 


V = / d 6 x 7 —> — i/3 / d 2 x_\_^/E = —i/3V E 


7 = V~9 -> ~IE = Vg = !>, 

J duj —f and uj —f uj n — n, 


(47) 


(48) 

(49) 


where V is defined first under eq. fl26|) . and the reason for u n = n is that, generally, u> n 
in the 5' 1 -compactification with a period R is given as oj n = 2 tt n/R, and in our situation 
R = 2tt. Correspondingly, the values £4 takes are discretized, and we change the notation 
of £4 as £4 —y fl n . 

As a result, the effective potential T defined as W E = /3 V E F with exp (iW) = 
exp(— W E ) can be obtained as 


r 


In the above, d 2 k± 


P 


7r u 


E 


' d 2 k± 

(M 


r-M 2 


dA 2 dXX sinh( 7 rA) 


K !xM 

A 2 -^ 


jdXX sinh( 7 rA) & J dk fcT. 


2nkdk (0 < k < oo), 


T = 



(50) 


(51) 


and we have evaluated the summation of n as 


E 


72— — OO 


1 


= f coth (it (/i — \/X 2 + /i 2 )) — coth (vr(/i + \/A 2 + /i 2 ))) 

yA 2 + /i 2V 2 

7T$. (52) 


Next, since we finally calculate the critical acceleration of the Bose-Einstein condensa¬ 
tion, we focus on the critical moment. The value of the chemical potential at the critical 
moment is given by m/a c , where a c is the critical acceleration. We explain this in what 
follows. 


12 










3.4 The chemical potential at the critical moment 


Using the effective potential T given in eq. (l5Ul) . the particle density d = — dT/dp can be 
written as 

d = -JL Jd\X sinh(vrA) ^ Jdkk^ + Q jdkk . (53) 

Here, to perform we show the /i-dependence of 4 by expanded it around p = 0 as 


with 


4 = 4 


94 

m=o 9/i 

m 2 

2 


fi =0 


ia 2 i< 


/i 2 + 


fi =0 


dA K iX (k, P ) - a K iX (np) 


'0 


A 2 =m 2 


/i + • • • 


94 
dp 
1 9 2 4 


_ dM 2 94 

n=o dp dM 2 


= 0 , 


2 dp 2 


/ 1=0 


fi=0 

1 9 / 9M 2 94 \ 

2 dp \ 9p 9M 2 J 


94 


= —a 


li=0 


dM 2 


A 2 =m 2 


Using the above expansion, we can write the particle density as 


d = — jdXX sinh(7rA) j"dkk ^ — 2a 2 /i &K 2 x (np) 


A 2 =m 2 


(54) 

(55a) 

(55b) 

(56) 


We can always confirm numerically that < 0 and 4 > 0. Further, we can see 4 > 0, 
since the integrand in 4 is given by the square and the integral direction is positive, as 
can be seen in eq. flbT]) . 

Hence, looking at eq. (l56]h we can say the following two things: When the chemical 
potential is increased with fixing all the other parameters, the particle density always 
grows up. On the other hand, when the acceleration is decreased with fixing all the other 
parameters, the particle density always decreased. 

Hence, when decreasing the acceleration from higher accelerations, where there is no 
condensation, say a = 0, with keeping the particle density to constant, the chemical 
potential should grow up to keep the particle density to constant. However, as can be 
seen from eq. (l22l) . there is the upper bound for the value the chemical potential can take, 
which is m/a. Hence, when the chemical potential reaches the upper bound in such a 
situation, a should start to grow up, where a plays the role of the condensation as written 
under eq. ffTTl) . 

Hence, we can see that the value of the chemical potential at the critical moment is 
m/a c . 


3.5 Calculation in the relativistic situation 

As mentioned under eq. (l52]h our analysis is the one at the critical moment. Hence, for 
the reason mentioned above, m/a c is assigned to the value of the chemical potential with 
a = 0 and p = a” 1 . 
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We put an assumption that m/a c is small. This means that the chemical potential in 
this study is also small. The situation that m/a c is small corresponds to the relativistic 
situation as mentioned in the introduction. 

In what follows, we use two symbols: /i, / x!. /i means the true chemical potential, where 
“true” means that it has been given at eq.((6]) and the derivative with regard to /i can act 
on. On the other hand, // means just a value m/a c . id' is not the chemical potential and 
the derivative cannot act on. 

Then, it turns out that we can expand <3? and T with regard to /i and /i' as 


$ = $ 0 + <h 2 /i 2 H-, 

T = T 0 /d ' 2 + T 2 ii 2 + ■ ■ • , 

where 

2 

$o = -r coth(7rA), 

A 

$2 = - 4 fcoth(7rA) + —T^y—-r(l - 27 tA coth(7rA)) V 

A^ V smh (7 tA) / 

T 0 = -T 2 = a 2 c K iX (k/a c ). 


(57a) 

(57b) 

(58a) 

(58b) 

(58c) 


At this time, substituting (I57aj) and (I57b|) into eq. ([30]) . T at the critical moment in the 
relativistic situation is given as 


r 


i 

2a c 7r 3 


dX A sin 



(59) 


4 The particle density and critical acceleration 

We now obtain the particle density d = —dY/djd using eq. fl59l) as 


d = 


-/ dXX sinh(7rA) dkk< h 0 T 2 

a c ^ Jo J 

2 a cf i r00 


7T J 

a^jd 

7T 2 


j dX cosh(7rA) J dkk K~ x (k/a c 
dX X coth (vrA), 


(60) 


where, in the above calculation, we have used sinh(7rA) To 1 

/ ttA 

dkkK 2 x (k/a) = — a 2 


= |cosh(7rA), and performed 
csch(7rA). 


Writing coth(7rA) = ^ + e2 ^i_ 1 in eq. (l60i) . we can see that the contribution from the 
1/2 diverges. A similar divergence appeared in a similar scalar field model with ours [IT], 
which was canceled by a mass renormalization. There also appeared the divergence in 
the Nambu-Jona-Lasinio model in the accelerating frame [15] [16], which was contained 
by a cutoff. We may be also able to show that the 1/2 can be canceled by some regu¬ 
larization. However there are some unclear points in its physical meaning. The source 
of the divergence can be thought clearly as the ultra divergence in the one-loop order. 


14 












Because our analysis can be regraded as an one-loop calculation from eq. (l24|) . and A can 
be regarded as some momentum from eq. flHFjh However, if we perform a mass renormal¬ 
ization by replacing m 2 with m q + dm 2 toward the m in the replacement of /i = m/a c to 
be performed right after this as the critical moment, it seems that the counter term 5m 2 
becomes independent of the acceleration. In this case, the counter term determined in 
the accelerating frame can survive even in the inertial frame, As a result, the consistency 
in the link between the accelerating and inertial frames in a model is unclear. On the 
other hand, if we perform the regularization by using a cutoff, how to determine the cutoff 
concretely is imclear for us at this moment. We do not examine this divergence problem 
in this paper, and handle it by simply ignoring the 1/2, considering that it is canceled by 
some regularization. 

As a result, by putting /i = m/a c as the critical moment, we can obtain the particle 
density as d = a 2 m/24:7r 2 . Hence we can obtain the critical acceleration as 



For the actual observation, we show the result in the MKS units as 



where d, and m in the last equation are the number of particle in a centimeter unit cube 
and the weight of particle measured by kilogram, respectively. 

5 Remark 

In this paper, based on the thermal excitation brought into by the Unruli effect, we 
have calculated the critical acceleration at which the Bose-Einstein condensation starts 
to occur, as in eq. flfiTT) . 

The system we have considered is an ideal gas composed of a complex scalar particle in 
the uniformly accelerating Minkowski space-time at zero-temperature. To this purpose, 
we have considered a free complex scalar field at finite density given in eq. (JTJ) in the 
Rindler space given in eq. (J4J) , and adopted an expression as in eqs. (H7a| and (11 7bh to 
describe the Bose-Einstein condensation. 

In our analysis, we have given the space-time integration in the action as the one 
without the p-integration to fix the acceleration and the Unruh temperature. This is 
different point from refs. [151 fTB] ITT] . If performing the analysis like them, there will be 
no difference in the final result. However, it is thought that we should describe explicitly 
that the space-time integration in the action is the one without the p-integration to make 
clear the physical meaning. 

To carry out the calculation, we have picked up the leading terms in the integrals 
§7D. As a result, we could obtain a result as in eq. (16 1 ]) . This manipulation means to 
assume that the complex scalar held considered in this paper has the critical temperature 
for the Bose-Einstein condensation in the temperature region fcgT/2 me 2 , which is the 
relativistic situation that the particle’s thermal motion energy is much higher than its 
static energy at the critical moment. 
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We can consider the analysis in another extremal region such as the non-relativistic 
situation in which the particle’s thermal motion energy at the critical moment is much less 
than the particle’s static energy as ksT /2 -C me 2 , where this T is around T c . However, 
the middle calculations appearing in the expansion in this limit are very complicated, and 
the analysis in this limit is technically difficult. We take this to a future work. 

Furthermore, the problem on the divergence in eq. dhHll has not been examined in this 
paper. We also take this to a future work. 

However, as mentioned in the introduction, there are papers concluding that the Unruli 
effect does not contribute to the restoration of the spontaneous symmetry breaking [18]. 
The system in these papers is not the finite density unlike us, and their conclusion would 
not be applied to our study readily. 

However, we can see intuitively that the situation which our study and the refs. m 
mm can lead is strange. Because, for a situation, it can be observed in different ways 
depending on the frame in which the observer stays. For example, for an observer in 
an inertial frame, it can be observed as a Bose-Einstein condensation state. But for an 
accelerating observer, it can be observed as an normal gas. Therefore they might be right 
after all. 

If we perform this study more, before that, it is thought that we should make clear 
whether their conclusion can hold or not in the finite density models like ours. 

In Ref. [24]. a critical temperature has been obtained in the same model with ours but 
the background space-time is taken to the Euclid space at finite temperature. Their result 
is ^/3d/m as written there, and we can see that there is difference from our result in the 
factor. We think that this discrepancy originates in the difference between the Rindler 
space that has the event-horizons and the Euclid space and what our calculation has been 
performed only in the right wedge of the Rindler space. 

We can see from our result that the critical acceleration rooted in the thermal excita¬ 
tion brought into by the Unruli effect has different behavior from the critical temperature 
obtained by the true thermal excitation in the flat space. We may expect from this that 
the observational data’s behaviors of the critical temperatures measured in the accelerat¬ 
ing and the flat static frames are different each other. It would be interesting, if we can 
find this discrepancy in the experiments. 

However, on the other hand, we think at last that the significance of this study would 
be what we could recognize a puzzle mentioned above. 
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